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Abstract 



Qh! We make a detailed analysis on the nucleon decay in the minimal super- 

symmetric 577(5) grand unified model. We find that a requirement of the 
unification of three gauge coupling constants leads to a constraint on a mass 
M Hc of color-triplet Higgs multiplet as 2 x 10 13 GeV < M Hc < 2 x 10 17 GeV, 
taking both weak- and GUT-scale threshold effects into account. Contrary 
to the results in the previous analyses, the present experimental limits on 
the nucleon decay turn out to be consistent with the SUSY particles lighter 
than 1 TeV even without a cancellation between matrix elements contributed 
from different generations, if one adopts a relatively large value of Mh c 
(> 2 x 10 16 GeV). We also show that the Yukawa coupling constant of color- 
triplet Higgs multiplet does not necessarily blow up below the gravitational 
scale (2.4 x 10 18 GeV) even with the largest possible value of M Hc . We point 
out that the no-scale model is still viable, though it is strongly constrained. 



1 Introduction 



The hierarchy problem has been the most serious problem in the grand unified theory 
(GUT) [|I|]. At present, the only feasible solution to this problem is to introduce the 
supersymmetry 0. Furthermore, the supersymmetric (SUSY) SU(5) model || is now 
strongly supported phenomenologically by the sin 2 9w measurement [[§] made at the LEP 
experiments ||. Once we regard the SUSY-GUT as a serious candidate of the physics 
beyond the standard model, a natural question is how we can test the model. The most 
striking consequence of the grand unification is the instability of nucleons. However, 
the nucleon decay via exchanges of X and Y gauge bosons is strongly suppressed as 
r~p oc Mqu T in the SUSY-GUT because of the large unification scale 

M GUT ~2x 10 16 GeV. (1.1) 

On the other hand, the nucleon decay via exchanges of color-triplet Higgs multiplet ||, 
which is suppressed only by MJ^ T , may still allow us to verify the model in the near 
future. 

The main purpose of this paper is to study the implication of the present experi- 
mental limits on the nucleon decay in the minimal SUSY SU(5) GUT (MSGUT). Sim- 
ilar analyses have been carried out by Ellis, Nanopoulos, and Rudaz ]7J, and later by 
Arnowitt, Chamseddine, and Nath || [| rather thoroughly. However, there has been no 
criterion given on how heavy the color-triplet Higgs multiplet can be. In this paper, we 
examine the experimental limits in the most conservative way, making the color-triplet 
Higgs multiplet as heavy as we can, allowed from the renormalization group (RG) analy- 
sis of the gauge coupling constant unification ||10|| . As a consequence, we find weaker con- 
straints than those given in the previous analyses. The authors of Refs. || have claimed 
that the data of nucleon-decay experiments at that time are already stringent enough 
so that the SUSY particles below 1 TeV are excluded unless there is a delicate can- 
cellation between the proton decay matrix elements from second- and third-generation 
contributions. On the contrary, we find that the present limits from the nucleon-decay 
experiments are still consistent with the SUSY particles below 1 TeV even without such 
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a cancellation. We also study a possible reach of the superKAMIOKANDE experiment. 
It will be shown that superKAMIOKANDE, together with LEP-II, is capable of covering 
most of the region with SUSY particles below 1 TeV, and hence it is highly expected to 
observe the nucleon decay at superKAMIOKANDE. It will be also stressed that more 
precise measurements on the gauge coupling constants, especially on that of QCD, will 
give a strong impact on the determination of the color-triplet Higgs mass Mh c . 

The paper is organized as follows. A brief review on the MSGUT is presented 
in Sect. 2, to summarize our conventions. We critically re-examine the analysis of 
Refs. |§ H] in Sect. 3. We find that the coefficients of the dimension-five operators 
are larger than theirs by a factor of 2. The decay rates for various modes are pre- 
sented. As pointed out in Ref. ||, there may occur a cancellation between second- and 
third-generation contributions. We present the partial lifetimes in terms of unknown 
parameters y tK or y t7T which represent the ratios of the third- to the second-generation 
contributions. In Sect. 4, we give an upper bound on the mass of color-triplet Higgs 
multiplet, requiring that the gauge coupling constants are unified. Then the present ex- 
perimental limits are examined in Sect. 5. There it is explicitly shown that the present 
data still allow for SUSY particles below 1 TeV, even without the cancellation between 
matrix elements mentioned above. The reach of the superKAMIOKANDE and the LEP- 
II experiments is discussed in Sect. 6. Sect. 7 is devoted to conclusions and discussions. 
An analysis on the dimension-six operators is presented in Appendix A. Appendix B 
summarizes discussions on the renormalization effects on the dimension-five operators. 
We improve the analysis given in Ref. J5J, but the difference turns out to be small. The 
chiral Lagrangian technique adopted to calculate the nucleon-decay matrix elements is 
described in Appendix C. 

2 Minimal SUSY 577(5) GUT 

In this section, we review the minimal SUSY SU(5) GUT (MSGUT) summariz- 
ing our conventions. We also clarify the origin of new CP-violating phases in Yukawa 
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coupling constants of color-triplet Higgs multiplet to matter multiplets. 

There are quite a few multiplets in the MSGUT. An adjoint Higgs multiplet £(24) 
breaks the SU (5) GUT group down to SU (3)c x SU (2) L x U(1) Y , and a pair of quintets 
77(5) and if (5*) contain doublet Higgs multiplets Hf, Hf in the minimal SUSY standard 
model as well as their color-triplet partners H c , He- The superpotential of this model 
is 



W = ^Tr£ 3 + -fVTiJ: 2 + A77 Q (£^ + W5 a p )H p 
3 2 

Uij 



(2-1) 



where the Latin indices i,j = 1,2,3 refer to families, and the Greek ones a,(3,^--- 
represent the 377(5) indices. The chiral superfields ip(10), 0(5*) are matter multiplets. 
Contents of the Higgs multiplets are 




-(3,2) 



(3*, 2) 




{He, H c , H c , Hj, Hf), 



H = {H c , H c , H c , H f , -H f), 



and those of the matter multiplets are 



* = 71 
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(2.2) 

(2-3) 
(2-4) 



(2-5) 



where all the matter multiplets are written in terms of the chiral (left-handed) super- 
fields. The chiral superfields u and d contain left-handed up-type and down-type quarks, 
u c and d c the charge conjugations of right-handed up-type and down- type quarks, e and 
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v left-handed charged leptons and neutrinos, and e c the charge conjugations of right 
handed charged-leptons. In the following, t Q = (u, d) and t L = (u, e) will denote chiral 
superfields of weak-doublet quarks and leptons, respectively. 

The SU (5) GUT symmetry is broken by a vacuum expectation value of the £ field, 



\ 



V 



giving masses to X and Y gauge bosons 



(2.6) 



M v = M x = M Y = h\/2g 5 V, 



(2.7) 



where is the unified 577(5) gauge coupling constant. The invariant mass parameter 
of H and H is fine-tuned to realize masslessness of Hj and Hf, while it keeps their 
color-triplet partners, H c and H c , superheavy as 



M, 



H C 



M- 



Hc 



h\V. 



{21 



The components Eg and E3 acquire the same mass 



M s = M S8 = .Uv, = -fV, 



(2.9) 



while the (physical) components S(3« j 2) and E( 3i2 ) form superheavy vector multiplets of 
mass My together with the gauge multiplets. The mass of the singlet component £24 is 
(l/2)fV. 

To analyze the dimension-five operators, we have to examine the Yukawa couplings 
of H and H to matter multiplets. An important question is how many independent 
parameters we have in the Yukawa couplings |]T2"fl . The Yukawa coupling constants h 1 ^ 
and f %3 in Eq. Q2.1| ) form a parameter space C 6 x C 9 , since h %3 is a symmetric matrix. 
The freedom of field re-definition is U(3) x {7(3), corresponding to the choice of the basis 
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of ipi and <pi. Thus the physical degrees of freedom of the Yukawa coupling constants is 
(6 + 9)x2-9x2 = 12 = 3 + 3 + 4 + 2. First two 3's stand for the eigenvalues for up- 
and down-type mass matrices, 4 for the Kobayashi-Maskawa matrix elements, and 2 for 
the additional phase degrees of freedom. We will parameterize the coupling matrices h % i 
and / u as 

h ij = tfe^Sv, (2.10) 

f' J = v; 3 p, (2.ii) 

with Vy being the Kobayashi-Maskawa matrix. Only two of the phases e lLpi are inde- 
pendent, and we can take 

(p u + <Pc + <Pt = 0. (2.12) 
In this parameterization, the corresponding bases of the matter multiplets are 

^ 3 *Qi = ( Ui , $ = (u h Vijdj), (2.13) 

ipi 3 e-^<, (2.14) 

ipi 3 V ifi % (2.15) 

0, 3 (2.16) 

<P % 3 % = (^, ei ), (2.17) 

in terms of the mass eigenstates u iy d^ u^, d°, z/j, e^, e\. Then the Yukawa couplings of 
Higgs to matter multiplets are given by 

W Y = tiQ i u c i H f + V* j f j QdpIf + fie c i Lji f (2.18) 

+^h i e i * i Q i Q i H c + VtftQiLjHc 
+h i V ij u c i e c j H c + e-^V^u-rdfHc. 

It should be clear from the above expression that the phases e t(fii would be completely 
irrelevant if He and He were absent. The phases appearing in the Yukawa couplings of 
He and H c cannot be absorbed by the field re-definition without affecting the couplings 
of Hf and Hf. As we will see later, these phases are important in the nucleon-decay 



amplitudes induced by the He and He exchanges. However, they are perfectly irrelevant 
to the nucleon decay caused by the X and Y gauge-boson exchanges. 

3 Dimension-Five Operators and Decay Rates 

In this section, we re-examine the previous analyses by Ellis, Nanopoulos, and Rudaz 
0j, and by Arnowitt, Chamseddine, and Nath || Several corrections to the formula 
and numerical factors are made, and as a consequence nucleon-decay amplitude turns 
out to be smaller than their result by a factor of 2. 

In the SUSY-GUT there are several baryon-number violating operators since it has 
many scalar bosons with color quantum numbers. Dimension-six operators induced 
by the X and Y gauge-boson exchanges are suppressed by 1/M GUT . These operators 
cause unacceptably large nucleon-decay rates in the minimal non-SUSY SU(5) GUT 



14] , but there is no problem in the MSGUT since Mqut is much larger than in the 
non-SUSY case. In fact, we have analyzed the nucleon decays caused by the dimension- 
six operators, and found that they are always suppressed compared to those caused 
by the dimension-five operators. A brief discussion on the dimension-six operators is 
presented in Appendix A. The dimension-five operators are much more dangerous. These 
operators are generated by H c and H c exchanges in the MSGUT, and is suppressed only 
by 1/Mgut P- The nucleon-decay amplitudes are obtained by dressing these operators 
by SUSY particle exchanges to convert scalar bosons to light fermions. Therefore, the 
nucleon-decay rates are sensitive to SUSY particle masses as well as the mass of He and 
He- It has been also noted that there may be baryon-number violating dimension-four 
operators [ 13| . However, they can be forbidden by imposing the i?-parity invariance, 



and we will not consider them in this paper. 

Let us now discuss the dimension-five operators that cause the nucleon decay. A 
supergraph is presented in Fig. 1. The operators can be written explicitly as 

W 5 = T^—Ve^VCtfiQiQJiQkLi) + ^h%e^V^f l (^ef)(u e k df) 1 (3.1) 



6 



where the contraction of the indices are understood as 

{QiQi){Qkh) = e aP Mdf '-dfvlXufa-dZvfi, (3.2) 
«$(i#J?) = e^v^epipfy. (3.3) 

with a, (3, 7 being color indices. Note that the total anti-symmetry in the color index 
requires that the operators are flavor non-diagonal (i ^ k). Therefore dominant decay 
modes in the MSGUT generally involve strangeness, like n, p — > Kv fl5| , 0j. 



The dimension-five operators will be converted to four-fermi operators at the SUSY 
breaking scale, by exchanges of gauginos or doublet Higgsino. However, the important 
contributions come only from the charged-wino dressing of the (QiQi)(QkLi) operators, 
and we will concentrate to this case. 

The exchanges of gluino, neutral gaugino and neutral Higgsino are small in general 
0, since they are flavor diagonal and hence suppressed by the Yukawa coupling constants 
of first and second generations appearing in the dimension-five operators. Though the 
gluino exchanges have stronger gauge coupling 0:3 than the wino exchanges, it will vanish 



completely in the limit where all the squark masses are degenerate fll5| . Since the high 
degeneracy is required to suppress the unwanted flavor changing neutral current |] the 
gluino exchanges turn out to be always small. The charged-Higgsino exchanges are 
also suppressed due to their small Yukawa coupling constants to the first or second 
generations.^ Thus the charged-wino exchanges give dominant contributions. Note that 
charged-wino dressing is impossible for the second operators {ulef){u%df) , since all the 
fields involved are right-handed fields. Though the left-right mixing due to the A-terms 
and superpotential \dW/dH\ 2 induces the wino dressing to the operators (u^) (u^df ) , 
their contribution is suppressed similarly to the charged-Higgsino exchanges since the 

*A high degeneracy is required to suppress flavor-changing neutral currents, especially between first 
and second generations. See Ref. |l6| . 

t There are contributions from the third generation to the charged-Higgsino exchange amplitudes. 
However, the bottom-quark Yukawa coupling constant f b is smaller than the SU (2) gauge coupling 



constant 52 unless tan/3// in Eq. (3.5) is extremely large. Therefore, the charged-Higgsino exchanges 



are most likely smaller than the charged-wino exchanges. 



7 



mixing is proportional to the Yukawa coupling constants. In the following we refer to 
the charged-wino simply as wino. 

By dressing the dimension-five operators with the wino exchanges, we will get four- 
fermi operators. The results depend on the masses of charginos, squarks, and sleptons in 
the loops. There is a mixing between wino and charged Higgsino, with the mass matrix 



M, 



chargino 



(3.4) 



V 



rriyj \/2mw cos (5h 

\f2m w sin 3 H /i 

Here, is a pure wino mass, /i a pure Higgsino mass, and Bh a vacuum angle of 
doublet Higgs scalars, defined by 

(H° f ) 

tan Bh = n • (3-5) 
( H f) 

The interaction of squarks and sleptons with wino is fixed by 

C = g2(u* L w + dL + d* L w~UL + v* L w + eL + e* L w~vi) + h.c, (3.6) 
giving the triangle diagram factor [f| 

g/K d) gi J — ) (^Tp) n ' (3 - 7) 

We have taken only the (1, 1) component of the chargino propagator, since the nucleon- 
decay amplitudes are dominated by the pure wino component in the chargino states. 
Though the integral Eq. ( |3.7| ) depends on the mass eigenvalues and the mixing angles, 
we have found that it is well approximated by the pure wino exchange, and then it can 
be given by 



fS { u,d)^f^^-J^^^4-^^^\. (3. 
2n 2n mf, -m^ - m~ m% mj - mi, mf 7 , 



This approximation can be easily justified if ^> mw, since then the off-diagonal 
elements in the mass matrix Eq. ( |3.4|) can be neglected. On the other hand, if ~ m-w, 
the off-diagonal elements cannot be neglected in general. However, if m^m^ 3> m^, 
the triangle diagram factor / in Eq. Q3.7D is simply given by 

(-^chargino) 11 TCl-fa , . 

J — 5 = 2~' (A^U 

mj, mj. 



and hence the above approximation is again justified. The region ~ ttlq ~ ~ 
/i ~ m w requires an exact treatment of the mixing, but this region turns out to be 
already excluded, and is irrelevant to our analyses. 

Notice that the nucleon decay rates depend on the SUSY particle masses only through 
the function /. It is useful to see the dependence of the function / on m^, nig, and 
m£. In the limit <C uIq ~ m^, / behaves as in Eq. ( [3.91 ). In the other limit 
rriu, ^> rriQ ~ mg, it behaves as 

1 777,?. 

ln^f. (3.10) 

rriw m z g 

These behaviors will be used to put bounds on these masses in section 5. 
The resulting four-fermi operators can be written down explicitly as 

x [(u?d!f)(d'j / v k )(f(u j , e k ) + f(u h d[)) + (c#X )(u]e k )(f(u h di) + /(dj, Vh)) 
+(d?u k )(dfu])(f( Ui , e k ) + f(u h d'j)) + (ufdf )(«7et)(/«, u 3 ) + /«, v k )) , 

(3.11) 

where the contraction of spinor indices are taken in each brackets (). Here we have 
assumed that the mixing between squarks is negligible. This is true in most of the su- 
pergravity models (for example, see Ref . [|l7jj ) which ensure the absence of flavor-changing 
neutral current. Notice that Eq. (|3.11|) is larger than that given in Refs. || || by a factor 



of 2. We suspect that this difference arises from an inconsistency of the normalization 
of their Yukawa coupling constants (See Eqs. (1.5), (1.6) in Ref. 0). 

We have to include three kinds of renormalization effects to perform quantitative 
analyses. First, the Yukawa coupling constants appearing in the Eq. ( |3.1|) are those eval- 
uated at the GUT-scale, and we have to calculate their magnitudes using the low-energy 
quark masses. Second, the dimension-five operators receive anomalous dimensions due 
to the wave-function renormalizations of the external lines. Third, the four-fermi oper- 
ators obtained after the gaugino- dressing will be further renormalized from the SUSY 
breaking scale down to 1 GeV. These three effects are first discussed by Ellis, Nanopou- 
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los, and Rudaz j7|. However, they dropped the SU(2) gauge interactions in estimating 
the first renormalization effects, which led to an overestimation by 50%. They also 
neglected the contributions from the top-quark Yukawa coupling. Since the top-quark 
mass is heavier than 90 GeV, these contributions, which appear in the wave-function 
renormalization of Higgs doublets, enhance the dimension-five operators by ~ 30%. 
However, corrections on these two points give roughly the same result as theirs. The 
details are explained in the Appendix B. 

After taking the renormalization effects into account, the nucleon-decay operators 
at 1 GeV are given by 



xe 



a/37 



[u?d'f)(d'?v k )(f( Uj , e k ) + f( Ui , <)) + (dfu?)(u]e k ) (f(u h d t ) + /(dj, v k )) 



+ (d' i a iy k )(dfu])(f(u t , e k ) + f(u h d')) + «df)(«7e fc )(/«, u 3 ) + /«, u k )) 



(3.12) 



where As(i,j,k) represents the short-range renormalization effect between GUT- and 
SUSY breaking scales depending on the flavor k, and Al the long-range renormal- 
ization effect between SUSY scale and 1 GeV. Here, the quark masses m u . and are 



defined at 1 GeV in the MS scheme 0. 

We first show the prediction of nucleon-decay rates for the dominant modes, n, p — > 
Kv^. The main contributions in the four-fermi operators Eq. ( |3.12j ) come from the terms 
(i = c, j = u, k = s) (proportional to m c m s ), and (i = t, j = u, k = s) (proportional to 
m t m s ). The relevant terms are given by 

C = P- rflit £ ^ ((^^ + ^){<Tv»)) 
M Hc m w sin zp H y ' 

x A s (c,u, s)m c e ilpc V cs V cd (f(c, fx) + /(c, d')) 
+A s (t,u,s)m t e i * t V ts V td (f(t, fi) + f(t, d'))] . (3.13) 

We have neglected the terms propotional to m u . Though the terms coming from the 
c-exchange can be computed precisely, the contribution of t-exchange is ambiguous due 
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to the unknown Kobayashi-Maskawa matrix elements for top quark fl8 |. In fact, the 
ratio of the t-contribution relative to the c-one [EL 



tK _ m t e^V ts V td A s (t,u,s)(f(t, fi) + f{t, d')) 
V ~ m c e i ^V cs V cd A s (c, u, s)(f(c, fi) + f(c, d')) 



ranges between 



0.096 < \y tK \ < 1.3, 



(3.14) 



(3.15) 



if we take the triangle diagram factors / to be common, and the top-quark mass m t to 
be 100 GeV.0 

Note that e lipt and e %LPc in Eq. ( 3.14j ) are independent of each other. We cannot 
measure these phases from the present-day experiments, since they are irrelevant to any 
obsevables as far as Hq and Hq are decoupled. Thus we are completely ignorant whether 
y tK is constructive or destructive to the c-exchange amplitude. We should regard this 
complex parameter free in the present analyses. 

If the modes n, p — > Kv^ have a cancellation between c- and t-exchange amplitudes, 
we have to study the other possible decay modes. Next-leading modes are n,p — > 
nu^, suppressed by the Cabbibo angle sin 2 Q c compared to the Kv^ modes. There are, 
similarly to the Kv^ modes, contributions from c- and i-exchange to these modes, and 
tliGir r£itjio 

u _ m t e^V t 2 d A s (t,u,s)(f(t, /x) + f{t, d')) 



y 



ranges as 



m c e^V^A s {c } u }S ){f(c, n) + f(c, d'))' 
0.041 < \y^\ < 1.7, 



(3.16) 
(3.17) 



with common / and m t = 100 GeV. Though y and y tv are correlated, we also regard 
y t7T as an independent parameter because of too large uncertainties in the Kobayashi- 
Maskawa matrix elements 

n.t-K T/ T/ 

(3.18) 







v td v cs 


ytK 




v ts v cd 



= 0.22 - 2.94. 

■^Note that m t in Eq. ( ^.14| ) is defined at the renormalization point 1 GeV. For example, m t 
270 GeV for m t = 100 GeV. 
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We find that a perfect "double" cancellation in Kv and iru modes is not possible, 
since y tK and y t7T have a non- vanishing relative phase coming from the CP- violating one 
in the Kobayashi-Maskawa matrix. Thus, we will not consider the "double" cancellation 
further in this paper, and take |1 + y t7T \ = 1 throughout. A "double" cancellation is 
logically possible only if the CP-violation is dominated by a non-KM mechanism. Even 
in the presence of such a "double" cancellation, there are decay modes which do not have 
such ambiguities (i.e., p — > K°fi + , r] fi + , 7r°/i + and n — > 7r~fi + ). The operators containing 
charged leptons have only a contribution from the up-quark Yukawa coupling constant. 
Therefore, the decay rates do not suffer from the ambiguity of a possible cancellation. 
However, the decay rates are very small since the up-quark Yukawa coupling constant 
is tiny. 

To obtain matrix elements at the hadronic level from the operators written in terms 
of the quarks fields, we adopt a chiral Lagrangian technique [19], p0~| . Details on this 
method is shown in appendix C. We present the results on the partial lifetimes of 
nucleons in Tables 1-3, in terms of the parameters 

/3,M Hc ,A s ,/3 H ,y tK ,y^, 

and the triangle factors /'s. A parameter (3 is the hadron matrix element parameter 
used in the chiral Lagrangian technique |H| ^2] 

Pu L (k) = e a ^(0\(d a y L )ul\p,k), (3.19) 

which ranges as 

P = (0.003 - 0.03) GeV 3 , (3.20) 

depending on the methods of the theoretical estimation. Due to the uncertainty of an 
order of magnitude, the predictions of nucleon partial lifetimes receive an ambiguity of 
two orders of magnitude. A more precise determination of f3 is strongly desired. 

Table 1 summarizes predictions on the nucleon partial lifetimes for the dominant 
modes n,p — > Kv. One sees that the partial lifetimes of neutron are a little shorter 
than ones of proton. This is because the former has a larger chiral Lagrangian factor. 
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Table 2 presents next-dominant modes n,p — > ttu. These modes are suppressed by the 
Cabbibo-angle sin 2 ^. Ratio of the decay rates into and v e is simply the squared 
ratio of strange- and down-quark masses. The decay rates into charged lepton fi + are 
listed in Table 3. These decay rates do not suffer from the ambiguity of a possible 
cancellation. However, the partial lifetimes are very long because of the small Yukawa 
coupling constant of up quark. 

The dimension-five operators given by Eq. ( |3.12j ) are larger than the expressions 
given in Refs. || Q . However, our final results shown in Tables 1-3 are smaller than the 
conclusion in Ref. by a factor of 2 which may originate in an inconsistency between 
their analytic formula Eq. (2.3) and its numerical evaluations Eq. (2.7) in Ref. |§. 

So far we have not discussed the decay rates of the modes containing v r . This is 
because we cannot make definite predictions for these decay modes since V* b has a large 
ambiguity. While the decay rates of v T modes can be as large as ones of P M mode if we 
take the largest possible V* b value, they can be also negligible with the smallest possible 
V* b value. In any case the decay rates can be only comparable to those into the v^, and 
hence the total decay rate into v is raised at most by a factor of two. This gives only 
a factor of s/2 stronger constraint on the squark masses, and we will not include the 
v T modes, hereafter. Once we know V* b more precisely, it is easy to incorporate the v T 
modes into the present analyses. 

Finally, we note that dimension-five operators depend sensitively on quark masses. 
We use the central values of current quark masses at 1 GeV which are estimated from 
the chiral perturbation theory and QCD sum rule in Ref. [p3fl . Quark masses of the first 
and second generations have large ambiguities. Especially, the strange-quark mass has 



a large error-bar (i.e., m s (l GeV) = 175 ± 55 MeV, given in the MS scheme). In our 
analyses we use m s (l GeV) = 175 MeV. 
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4 Constraints on the GUT-scale Mass Spectrum 



Since the nucleon-decay rates due to the dimension-five operators are proportional to 
the inverse square of the color-triplet Higgs multiplet mass M^,, it is very important to 
determine it by some means. In the previous analyses J7L EL M , the authors have chosen 
Mh c = (1 — 2) x 10 16 GeV ad hoc. However, we have shown recently flOf , that one 



can obtain limits on the GUT-scale mass spectrum in the MSGUT, just by requiring 
the unification of three gauge coupling constants. In particular, we have derived the 
upper bound on Mh c without any theoretical prejudice. A theoretical requirement 
that the Yukawa couplings remain perturbative below the gravitational scale Mp / y/8n 
(2.4 x 10 18 GeV), poses further constraint on the GUT-scale mass spectrum. 

We first discuss the renormalization-group (RG) evolution of three gauge coupling 
constants. It was shown in Refs. E3, E31 that the simple step-function approximation is 



accurate for supersymmetric theories, justified in the "supersymmetric regularization" 
DR scheme . To illustrate how the GUT-scale spectrum receives constraints, we first 
discuss the one-loop RG equations. After that we include the two-loop corrections. 

The running of three gauge coupling constants in the MSGUT can be obtained easily 
at the one-loop level as 

a - Vz) = a -. (A) + _L{(_ 2 _^) ln !^ + ( _ 9 + 2JV9)ln A 

A , A n A 1 , , 

-4 In— - + 31n — +m— — \, 4.1 

\m z ) = a -(A) + f {(-|-|AT ff -|)ln!^ + (-6 + 2Ar fl + l)l n A 

27T L V 3 3 6/ m z m z 

-6 In— - + 21n — }, (4.2) 



a 2 



M v My 



'V 

a - Amz) = ^(A) + ^{O s -i)ln^f + (27V 9 + ?).„A 

Here, the scale A is supposed to be larger than any of the GUT-scale masses. The 
number of generations N g is three, and we have assumed a common mass msusY for all 
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the SUSY particles and for the scalar component of one of the Higgs doublets. A mass 
of the other doublet Higgs boson is taken at mz- By eliminating a^ 1 from the above 
equations, we obtain simple relations: 

(3a--2a 3 - 1 -ar 1 )(-z) = i In ^ - 2 In , (4.4) 

Z7T I 5 mz mz J 

fr -1 o -1 o — 1\ / x 1 L, M£M s m 5C/5 yl 

5a, — 3a 2 —2a, ](mz = — < 12 In — —s horn >. (4.5) 

2n [ m z mz J 



The Eqs. (|4.4|) , ( |4.5| ) imply that we can probe the GUT-scale mass spectrum from the 
weak-scale parameters (i.e., gauge coupling constants and mass spectrum of the SUSY 
particles). § Especially, Mh g is determined independently of My and Ms. Eq. ( |4.5|) 
determines a combination of the vector and adjoint-Higgs masses (MyMs) 1//3 , and we 
will call it as "GUT-scale" M GUT = (M^) 1 / 3 , hereafter.Q 

So far we have assumed a common mass msusY for the SUSY particles, but the mass 
splitting among the SUSY particles is also important to determine the GUT-scale mass 
spectrum. To avoid unnecessary complications, we restrict ourselves to the minimal 
supergravity model |T7|1 , where the SUSY-breaking mass parameters at the weak-scale 
can be determined from a small number of parameters at the Planck scale, by using the 



RG equations [ 28[ . Therefore, the squark and the slepton masses are given by 



ml = m 2 + 6.28M 2 + 0.35m| cos 2/3#, 

m l = m 2 + 6.28M 2 - 0.42m| cos 2(3 H , 

m\ c = m 2 + 5.87M 2 + 0.16m| cos 2/3 H , 

m\ = m 2 + 5.82M 2 -0.08m|cos2/? H , (4.6) 

m \ = m 2 + 0.52M 2 - 0.27m| cos 2(3 H , 



Ht was claimed in Ref. that the threshold corrections at the GUT-scale is so large that one cannot 
predict the SUSY-breaking scale even if measurements on 013 become much more precise. This high 
sensitivity on GUT-scale mass spectrum implies that one can probe it through precision measurements 
on the weak-scale parameters. Therefore, our result is consistent with their claim. 

^This "GUT-scale" Mqut does not necessarily correspond to the scale where all three gauge coupling 
constants meet. 
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mj = m 2 + 0.52M 2 + 0.50m| cos 2/3# , 
m 2 c = m 2 + 0.15M 2 -0.23m|cos2/3 H , 

in terms of the universal scalar mass m and the gaugino mass M at the GUT-scale. We 
have neglected the contributions from the Yukawa couplings to the renormalization of 
the particle masses. Also, the gaugino masses at the weak-scale are given by 



a 5 (M GUT ) 

m. = -f^M, (4.7) 
a 5 {M GUT ) 

a 3 (m z ) 
a 5 (M G[/T ) M ' 



where and represent masses of bino and gluino, respectively. 

The effect of the mass splitting can be taken into account by replacing \rt(msusY l m z) 



in Eqs. (gj) , (gj ) as 



in Eq. (OH, and 



21n^^ — > 4 In ^ + ^ In 2 ~ ~ ln ~ ~ ~ ln ~ 

rriz rriyj 5 m%m~ 5 5 



81n m^_ 41n m, +41n m, + (4 . 9 ) 

in Eq. (|4.5|) . Two doublet Higgs bosons are assumed to have masses at m# and mz, 
respectively. The symbol m~ h represents a mass of doublet Higgsino. We have neglected 
the mixings among gauginos and doublet Higgsino. 

For the time being we will restrict ourselves to the case where the universal scalar 
mass dominates the SUSY breaking {i.e., m ^> M). The terms ln(m|cm| c mgc / 'm^rri^) 
in Eq. ( J4.8|) and In (mj~/mecm„c) in Eq. ( |4.9| ) are negligibly small. The term \n(mg/mw) 
stays constant, since mg/m^ = a 3 /a 2 ~ 3.5. The dependence on m H in Eq. (|4.8|) is 
weak due to its small coefficient, and we set m# = 1 TeV. Therefore, we find that Mh c 
depends mainly on the Higgsino mass m~ h , and Mqut on the product of gaugino masses 
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mgiriyj. We have also examined the constraint on Mh c in the no-scale model |29j {i.e., 
m = 0), and found that the difference is negligible. 

Now we are at the stage to derive the GUT-scale mass spectrum from the above RG 
analysis. In our numerical calculation, we use the one-loop RG equations for the weak- 
and the GUT-scale thresholds, and the two-loop ones between these two distant scales. 
The two- loop RG equations in the minimal SUSY standard model are [24| 



1 



16n 



where i, j = 1, 2, 3, and 
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(4.11) 



N. 



H, 



(4.12) 



Here, N Hf is the number of doublet Higgs multiplets (Njj. = 2). The threshold cor- 
rections at the two-loop level are expected to be small, since their mass splittings only 
within the same order of magnitude do not produce large logarithms. As the input 
parameters, we use the MS gauge coupling constants at the Z-pole given in Ref. |30 . 
a = 127.9 ± 0.2, sin 2 B w = 0.2326 ± 0.0008, and a 3 = 0.118 ± 0.007. However, the 



use of the simple step-function approximation is only justified in the DR-scheme. Since 
we employ the simple step-function approximation, we have to convert these coupling 



constants at the Z-pole into the DR-scheme by 

1 1 Ci 



a 



DR 



a 



MS 



12tt' 



(4.13) 
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where d = 0, C 2 = 2, and C 3 = 3 |I . 



Combining all the above discussions, we find that M Hc is constrained to the range 

2.2 x 10 13 GeV < M Hc < 2.3 x 10 17 GeV, (4.14) 

and the "GUT-scale" is tightly constrained asjJJ] 

0.95 x 10 16 GeV < (M^M E ) 1/3 < 3.3 x 10 16 GeV, (4.15) 

for 100 GeV < rrig < 1 TeV. The allowed range of Mjj c is much less constrained, 
because of the small gauge-group representation for the Higgs multiplets. The large 
ambiguity of Mh c comes mainly from the uncertainty in the strong coupling constant a^, 
and the prediction on the nucleon decay will be drastically improved if the uncertainty 
diminishes. In Fig. 2 we present the GUT-scale spectrum derived from the present gauge 
coupling constants, and also that expected if the error-bar of 0*3 is reduced by a factor 
of 2 with the same central value. The importance of more precise measurements on 0:3 
should be clear from the figure. 

When one uses these constraints, Eq. (fO~4]) and Eq. Q4.15 ), one needs to pay two at- 



tentions. First, we have taken only one standard deviation for gauge coupling constants. 
If we allow two standard deviations, allowed region of Mh c spreads to both ends by two 
orders of magnitude, loosing any practical limits. On the other hand, the "GUT-scale" 
is still constrained tightly. Therefore, we restrict our analyses to only one standard 
deviation as in Eqs. ( [4.14|) , (|4.15|) .f"l Second, we have used the RG equations at the 



two- loop level. One may be concerned for whether three- loop corrections are important. 
The difference in Mh c between the one- loop and the two-loop results is a factor of 30, 
which is a very small factor compared to the large ratio Mqut/^z appearing in the 
II The bound on M GUT quoted in Ref. (TJ, 0.90 x 10 16 GeV < M GUT < 3.1 x 10 16 GeV, includes a 



minor mistake, and should be replaced by Eq. (4.15) 



**In other words, it is still possible to raise Mh c up to near the gravitational scale if we allow two 
standard deviations. However, we tentatively take this one-sigma bound seriously to perform further 
analyses. It should be noted that nucleon decay can be generated at observable rates for reasonable 
range of parameters, even with the color-triplet Higgs of mass at the gravitational scale. 
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solutions of the RG equations. Thus we expect that the three-loop corrections are much 
less than 0(1). 

Although we have concentrated to the RG analysis on the gauge coupling constants 
to determine the GUT-scale mass spectrum, we will obtain further constraint on the 
mass spectrum from the RG analysis on Yukawa coupling constants. As shown in 
Eq. (|2.8| ), M Hc is given by a Yukawa coupling constant A between H c , H c and E, 
which is not known. On the other hand, the mass My is determined by the SU(5) 
gauge coupling constant g§, whose strength is known by the RG analysis. A large 
mass splitting My <C M Hc requires that the A is very large compared to g 5 . Thus 
the applicability of the perturbation theory restricts the mass splitting to be not large. 
The same argument applies to the mass Ms, which originates in a self-coupling of the 
adjoint-Higgs as seen in Eq. ( |2.9|) . 

A constraint arises by requiring that those Yukawa coupling constants do not blow- 
up below the gravitational scale, Mp/\/8ii = 2.4 x 10 18 GeV. The running of the Yukawa 
coupling constants in Eq. (|2~T| ) are described by the RG equations, 



where h f is the Yukawa coupling constant to top quark. The conservative limit on A 
can be obtained in the case / = h f = 0. A numerical study shows that the mass Mh c 
is limited from above, 




(4.16) 



(4.17) 



(4.18) 



(4.19) 



M Hc = ^^M v < 2.0M V . 
A similar limit on Ms can be obtained with A = h l = 0, 



(4.20) 



Ms 



/ 



My < l.SMy. 



(4.21) 



2V2g 5 
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One may feel uneasy about the assumption that there is no new physics between 
the GUT-scale and the gravitational scale. We have examined the above analysis again 
requiring that the Yukawa coupling constants do not blow-up below 10 17 GeV. This 
requirement relaxes the constraint Eqs. ( 4.20Q , ( [4.21 ) at most by a factor of 2. Therefore, 
in the following calculation we use Eqs. (|4.20j) , (|4.21| ). 

We obtain constraints on My and Ms separately, combining the discussions above. 
The limits Eq. J05|) and Eq. flOID give 



M v > 0.78 x 10 16 GeV 



M s < 4.9 x 10 16 GeV. 



(4.22) 
(4.23) 



Eq. (|4.22 ) will be used to put limits on the dimension-six operators in the following 
sections. 



5 Present Limits on Dimension-Five Operators 

In this section we combine the predictions obtained in the previous sections with the 
results of the nucleon-decay experiments to see the present status of the MSGUT. We 
find that the present lower bounds on the nucleon partial lifetimes are still consistent 
with the SUSY particles below 1 TeV if one adopts a relatively large value of M# c 
(> 2 x 10 16 GeV). In Table 4, we have listed the experimental lower limits on the 
partial lifetimes of nucleon |18[ . 



The most dominant decay mode by dimension-five operators is n — > z/ M , as shown 
in Table 1. This mode dominates slightly over the similar decay mode, p —> K + u^, 
because of the chiral Lagrangian factor. These decay modes, however, have an am- 
biguity in the parameter y , which is the relative ratio of the t-exchange to the c- 
exchange contributions. Furthermore, the parameter y tK contains an unknown phase 
factor e 1 ^ - ^, which cannot be determined from any low-energy experiments. In fact, 
Arnowitt, Chamseddine, and Nath || have shown a possible cancellation between the 
second- and third-generation contributions. However, if the combination |1 + y tK \ de- 
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creases, the modes n,p — > Kv^ cease to be dominant. Then the experimental limit on 
the other modes n,p — > 7rz/ M become important. This interchange occurs at 

\l + y tK \ = 0.40, (5.1) 

when |1 + y t7T \ = 1. The parameters y tK and y t7T are correlated as clear from their 
definitions (Eqs. (|3.14j) , (|3.16| )). However, we find it impossible that |1 + y tK \ and 
|1 + y t7r | are canceled out simultaneously, as explained in section 4. We take |1 + y tn \ = 1 
throughout. 

In Fig. 3, we show the lower bound on Mh c derived from the present nucleon-decay 
experiments by varying |1 + y tK \. We choose other parameters such that the nucleon 
lifetimes become as long as possible {i.e., ttLq = = 1 TeV, = 45 GeV, and 
tan/3# = 1). In this figure, the upper horizontal line corresponds to the maximum value 
[Mu c = 2-3 x 10 17 GeV) in Eq. ( 4.14 ). There are two curves representing the lower 



limit on Mh c obtained from the experimental limits on the nucleon lifetimes. The upper 
curve corresponds to the case of the hadron matrix element (3 = 0.03 GeV 3 , and the 
lower curve to the case of (3 = 0.003 GeV 3 . The smaller hadron matrix element (3 gives 
weaker constraint as expected. Thus, the conservative lower bound on Mh c from the 
nucleon-decay experiments is 

M Hc > 5.3 x 10 15 GeV. (5.2) 

We illustrate how the lower bounds on M# c depend on the SUSY breaking param- 
eters, the wino mass and the squark mass itiq in Fig. 4, assuming ~ T^gl^l I n 
this figure the dashed line shows the dependence on when we choose the most con- 
servative set of parameters, ttiq = mi = 1 TeV, tan/^ = 1, |1 +y tK \ < 0.4, As = 0.67, 
and (3 = 0.003 GeV 3 . The lower bound on M Hc rises linearly on m^, in the region 
rriyj < 1 TeV. However, the lower bound decreases as m^ 1 beyond 1 TeV, and we do 
not obtain an upper bound on with this conservative choice of parameters. The 

^The situation does not change even if we allow mass splittings between sleptons and squarks (say, 
TO£ <C tUq), since the denominator in Eq. ( |3. 9| ) is dominated by the heavier mass, tUq. 
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dash-dotted line shows the dependence on ttiq again for the most conservative case, 
= 45 GeV, tan/3 H = 1, |1 + y tK \ < 0.4, A s = 0.67, and (3 = 0.003 GeV 3 . The lower 
bound on M Hc goes down as l/ m ^> leading to the lower bound on uiq, 

niQ > 150 GeV. (5.3) 

We also show the dependence on tan/3# in Fig. 5. We have fixed ttIq = mi = 1 TeV, 
= 45 GeV, (3 = 0.003 GeV 3 , and \l + y tK \ = 1.0 or \l+y tK \ < 0.4. The lower bound 
on Mh c is proportional to 

1 1/1 



+ tanP H . (5.4) 



sin 2 (5 h 2 ytan/?H 

When tan/3# 3> 1, the dependence is almost linear. We find a constraint on tan/?#, 

tan/?^<85, (5.5) 

which is, however, much weaker than tan (3h < 40 obtained from the requirement that 
the Yukawa coupling constant f b for bottom quark remains in the perturbative regime 
below the "GUT-scale" .0 

Taking Mh c as heavy as possible given in Eq. ( |4.14| ), we obtain limits on the masses 



rriu, and mQ. The allowed region is shown in Fig. 6. Here we have taken |1 + y \ = 1 
and mi ~ mQ. The present experimental limits on the wino and the squark masses from 
direct-search experiments at LEP |32[ and CDF |3^] are shown for comparison. Since 
the decay rate behaves like (m^/m^) 2 in the region m^ m^, the lower bound on m^ 
behaves like m^/ 2 . In the other extreme, m^ ^> m^, the decay rate goes like l/m%, 
and around m^, ~ 10 5 GeV the constraint on m^ from the nucleon-decay experiments 
becomes weaker than that from the CDF experiments. We see that the "natural" mass 
region < 1 TeV for the SUSY particles still survives the nucleon-decay experiments. 

Though the authors of Ref. claimed that the present limits on the nucleon decay 
are stringent enough to exclude the SUSY particles lighter than 1 TeV in the absence 

^The tan (3h has to be larger than 0.5, since otherwise the top-quark Yukawa coupling constant will 
blow up below the "GUT-scale" with m t > 90 GeV. 
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of a delicate cancellation between matrix elements of the dimension-five operators (i.e., 
(3 ~ 0.003 GeV 3 , |1 + y tK \ ~ 0.2), we see now that there is a wide allowed range. This 
is mainly because we use the maximum value of Mh c (2.3 x 10 17 GeV) given from the 
RG analysis while they chose M Hc — M GUT just by hand (i.e., 2.0 x 10 16 GeV). 

We show similar limits from the ttv^ mode in Fig. 7. As discussed above, these decay 
modes become dominant if |l + ?/ tj:s: | < 0.4. The limit is weaker compared to the previous 
case without the cancellation as shown in Fig. 6. 



We have taken the largest value of Mh c in Eq. ( f4.14|) in Figs. 6,7. This value requires 
My larger than 1 x 10 17 GeV in order to satisfy the requirement Eq. ( |4.20| ), which leads 
to M s smaller than 3 x 10 15 GeV. Though this case needs a mass splitting among the 
heavy particles, it is still within two orders of magnitude, and it is completely acceptable 
phenomenologically. 

In the minimal supergravity model, the SUSY particle masses are determined mainly 
by the universal scalar mass and the gaugino mass. In models where the universal 
scalar mass dominates the SUSY breaking parameters, squarks and sleptons are almost 
degenerate, that are larger than wino. Therefore, models of this type are preferred. In 
the no-scale model the SUSY breaking is dominated by the gaugino mass |2{J, which 
results in a definite prediction of the SUSY particle masses (nig ~ rrig ~ 3m^ ~ 3m^) 
f35fl . Since all the mass parameters in the triangle factor / defined in Eq. ( |3.8|) are 
proportional to m^, f behaves as . This enables us to derive the lower bound on 
from the nucleon-decay experiments. We have found that lower limit on is 70 GeV 
(equivalently, mQ > 210 GeV), in the most conservative case (i.e., (3 = 0.003 GeV 3 , 
M Hc = 2.3 x 10 17 GeV, tan(3 H = 1, and |1 + y tK \ = 0.4) .0 Thus, we conclude that the 
no-scale model is still surviving.^] 

Finally, we briefly discuss the nucleon decay caused by the dimension-six operators, 
which come from the X and Y gauge-boson exchanges. We show the decay rates by the 

*In the absence of the cancellation (i.e., |1 + y \ = 1), the corresponding limit is > 180 GeV. 
tThe wino mass in this region is consistent with the limit < 300 GeV |55| in the radiative 
breaking scenario of the electroweak symmetry in the no-scale model. 
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dimension-six operators in Appendix A. We find that dimension five operators always 
dominate dimension-six operators, in agreement with the naive expectations. The ratio 
(R T ) of the partial lifetime of the decay n — > K°i?^ via the dimension-five operators to 
that of the decay p — ► 7r°e + via the dimension-six operators is 



where all parameters besides Mh c were chosen so that R T becomes as large as possible. 
The ratio of Mu c to My is smaller than 2.0 from Eq. (|4.20p , and R T is always smaller 
than ^. The MSGUT predicts that the dimension-five operators should be observed 
earlier than dimension-six operators. 



6 Future Tests on the Minimal SUSY 577(5) GUT 



Now we examine how stringent constraint on MSGUT we can obtain from the nucleon- 
decay experiments in the near future. The superKAMIOKANDE experiment will push 
up the lower bound of the nucleon lifetime by a factor of 30. Meanwhile, the LEP-II 
experiment will be able to find wino below mz- Since the larger means faster nucleon 
decay as we discussed in previous sections, the combination of these two experiments 
can put more stringent limits on the dimension-five operators. 

In Figs. 8 and 9, we show the expected limits on the dimension-five operators if 
superKAMIOKANDE does not observe the Kv^ and ttu^ decay modes, and if LEP-II 
does not discover wino below mz- For the Kv^ decay modes, superKAMIOKANDE and 
LEP-II will be able to exclude most of the region with the SUSY particles below 1 TeV, 
even with the smallest hadron matrix element, (3 = 0.003 GeV 3 . In the case where 
the 7rz/^ mode is dominant, the region below 1 TeV is almost closed leaving a little 
window. Thus one can see that the LEP-II and the superKAMIOKANDE experiments 
are extremely important for testing the MSGUT. 



R T = 




(5.6) 
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One may be concerned for a little region below 1 TeV which may be left by these 
experiments. This window is open because of the large maximum value of Mh c obtained 
from the gauge coupling unification. If the error-bar of is reduced by a factor of 2 with 
the same central value in the future, the maximum value of Mh c becomes 6.1 x 10 16 GeV, 
and one will be able to close the window completely. In Fig. 10 we demonstrate the case 
where the error-bar of is reduced by a factor of 2. If one wishes to test the MSGUT 
completely, it is quite important to reduce the error-bar of 0*3. 

At the end of this section, we see whether the nucleon decay via the dimension- 
six operators can be observed in the near future. From Eq. (|4.22|) and Eq. ( |A.2 ), the 
theoretical lower limit of the partial lifetime of the decay p — > 7r°e + is obtained as 

r(p -> 7r°e+) > 4.1 x 10 33 ( °' Q3 ^ eV3 ) years, (6.1) 

in terms of the hadron matrix element a (=0.003-0.03 GeV 3 , see Appendix A). Since 
superKAMIOKANDE is expected to reach up to r{p — ► 7r°e + ) ~ 10 34 years f36fl , there 
is a possibility to observe the nucleon decay via the dimension-six operators in the 
MSGUT. 



7 Conclusions and Discussions 

We have analyzed the nucleon decay in the minimal SUSY 577(5) GUT (MSGUT) in 
details. First, we have studied the GUT-scale particle spectrum using the RG analysis, 
and found a maximum value of Mh c to be 2.3 x 10 17 GeV. Then, we have studied 
the nucleon partial lifetimes with the largest possible Mh c - We have found that the 
present nucleon-decay experiments are still consistent with the MSGUT, even without 
a cancellation between the matrix elements from exchanges of squarks in different gen- 
erations. We have emphasized the important role of precise measurement of the gauge 
coupling constants, especially the QCD coupling constant 03 to determine the mass 
Mh c ■ We have also stressed that the combined information of the lower bound on the 
chargino mass (LEP-II) and on the nucleon lifetimes (superKAMIOKANDE), will give 
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a stringent constraint on the MSGUT. 

It deserves to mention that the nucleon decay prefers supergravity models where 
the SUSY breaking parameters are dominated by scalar masses, rather than by gaugino 
masses. For example, the no-scale model, where the SUSY breaking parameters come 
only from the gaugino masses, is strongly constrained from the nucleon-decay experi- 
ments, though it is still viable. Increasing the limits on the gaugino masses will have 
strong impact on the MSGUT phenomenology. 

It is important to see whether predictions on the dimension-five operators become 
drastically altered by modifying the MSGUT. In the MSGUT, there is a prediction of 
mass relation, m& = m T , m s = m^, and uid = m e , at the GUT-scale. It is known 
that, though mb = m r is consistent with observations if top quark is not too heavy ||37|| , 
m s — n^n and vn^ = m e are not. Therefore, the modification of the Yukawa coupling 
structure is needed. Georgi and Jarlskog proposed that mj = m T , 3m s = m M , and 
rrid = |m e at the GUT-scale, introducing a 45 dimensional Higgs scalar f38fl . Also, 
Kim and Ozer proposed to make an "effective" 45 dimensional Higgs scalar by using 
higher dimension operators j^Jj. These modifications may produce more uncertainties 
in the Yukawa couplings of Hq and Hq to matter multiplets. However, we have checked 
that these models receive at most a few times stronger constraint, and hence the main 
conclusion in the present analyses does not change qualitatively. 

It has been argued that the introduction of a Peccei-Quinn symmetry may be able 
to eliminate the dimension- five operators |6j . However, the requirement of the coupling 
constant unification does not allow us to introduce new multiplets of arbitrary masses. 



We have shown in a separate paper |4(| that a suppression of the dimension-five operators 
by introducing the Peccei-Quinn symmetry cannot be stronger than in the MSGUT. 

Another interesting result of our RG analysis on the GUT-scale spectrum is that one 
can raise the grand unification scale up to the gravitational scale, by lowering the Ms 
down to the intermediate scale (10 11-12 GeV). In this case the heavy vector multiplet as 
well as He and He lie at the gravitational scale, if one allow for two standard deviations 
in Eq. ( |4.14| ). These light Eg and E3 require an extremely small Yukawa coupling /E 3 , 



26 



leading to a very flat potential of E. The possibility that the three gauge coupling 
constants meet at the gravitational scale has come out with a seemingly accidental 
parameter tuning. However, this may suggest a completely different underlying physics 



(non-GUT) like the superstring theory [41 
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Note added 

After completing this work, we have received preprints by R. Arnowitt and P. Nath 
In these works, they also derive a constraint on M# c , requiring the Yukawa coupling 
constants to remain in the perturbative regime. The difference between their constraint 
Mh c < 3 My and our Eq. ( 4.20Q is due to different requirements. Namely, they impose 



that the Yukawa coupling constant A should not blow up below 2M# C , while we impose 
that below the gravitational scale. However, they are not aware of the possibility to 
raise My by lowering Ms, and hence obtain a smaller upper bound on M# c . This leads 
to the opposite conclusion on the no-scale model, where they claim it to be excluded 
while we have found it still viable. 



Appendix A Decay Rates via Dimension-six Oper- 
ators 

In this appendix we analyze the nucleon-decay rates via dimension-six operators. The 
dimension-six operators are caused by exchanges of X and Y gauge-bosons or color- 
triplet Higgs scalars. The color-triplet Higgs scalars interact with matter only by small 
Yukawa coupling constants, and also its mass should be larger than 5 x 10 15 GeV (see 
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Eq. (|5.2|) ). Thus the dimension-six operators induced by the color-triplet Higgs scalar 
exchanges are negligible. 

The dimension-six operators via the X and Y gauge-boson exchanges are dominated 
by generation-diagonal decay modes, and the decay modes containing strangeness are 
suppressed by the Cabibbo angle sin 9c- We discuss here only the dominant decay mode 
p — > 7r°e + . The amplitude by the X and Y gauge-boson exchanges is 

C x ,y = - A « e< *"]^" c «^r + (1 + iK.il^CdgtiS)^^)) . (A.l) 



The renormalization factor Ar calculated by the authors in Ref. |43| is Ar = 3.6. 
Since the decay rate of this mode caused by the dimension-five operator is extremely 
suppressed, the observation of this decay mode would suggest the presence of the X and 
Y gauge-bosons exchanges. We calculate the decay rate using the chiral Lagrangian 
technique, as 



T 



M v \ 4 / 0.003 GeV 3x 2 



( P A+) = 11 x 10- x ( w ^ j ) \-^—\ y-s. (A.2) 



Here, a is the hadron matrix element 

au L (k) = ea^y(0|(da«fl)«lb, k), (A. 3) 



whose absolute value is the same as (3 (i.e., \a\ = \(3\) [^T|]. It is straightforward to 
compare it with the nucleon-decay rates via the dimension-five operators. The result is 
given in Eq. (|5.6|). 



Appendix B Renormalization Factors 

In this appendix, we present several formulae necessary to compute the renormalization 
factors which have appeared in the section 3. There are three kinds of renormalization 
effects. First, we have to derive the Yukawa coupling constants of He and He from the 
observed quark masses. Second, the dimension-five operators derived at the GUT-scale 
receive anomalous dimensions from the wave-function renormalizations of the external 



28 



fields. Third, the four-fermi operators dressed at the SUSY breaking scale should be 
renormalized down to 1 GeV. In the text the short-range renormalization factor between 
the SUSY breaking and the GUT-scales is denoted as As, and the long-range renormal- 
ization factor between 1 GeV and the SUSY breaking scale as A^. Though As and Al 
are estimated by the authors in Ref. 0j, we need minor corrections to their calculation 
of As- We demonstrate the derivation of As, and also comment on A^. 

Since the Yukawa coupling constants in Eq. ( |3.1| ) are those given at the GUT-scale, 
we have to calculate the values of the Yukawa coupling constants by solving the RG 
equations from the SUSY breaking scale up to the GUT-scale. This was done in the 
Ref. 0. Since the Yukawa couplings Hf and Hf are F-terms, what we have to compute 
is only the wave-function renormalizations of each chiral superfields thanks to the non- 
renormalization theorem of the F-terms [34]. The authors of Ref. [Q have given the 
following general formula for the running of arbitrary F-terms. Any F-terms can be 
re-written in the following form, 

W = ^h^kMm, (B.l) 

where 0's are chiral superfields, and the coupling constant h klm are supposed to be 
completely symmetric under the interchange of the indices k, I, m. If the fields </>fc's 
belong to the representation R k of certain gauge group i, the coupling constants h klm 
follow the RG equation 

dm i 

6 k k ,h k ' lm + 6\,h kVm + 9Zh klm '] , (B.2) 



dh klm 



dln/x (47r) 2 
where the constants #'s are given by 

Q k , = -2Cl(R k )g*5 k , + \h k ™h\ lvq . (B.3) 

The gauge coupling constant is denoted by gi, and the second Casimir of the represen- 
tation R k by C^iR 1 *)- If the coupling constant h klm can be neglected in the expression 
of 8, the RG equation can be easily integrated to the form 

/ ( \\Ci{R r )/bi 

h kim (fi) = h kim (/j ) n n Ptt • (B-4) 

r=k,l,m i=l,2,3 V a AH I 



29 



Here, the coefficient of the /3-function b{ is defined as 

da^ 1 bi 
d In jj, 2rr 



(B.5) 



The explicit forms of frj's in the minimal SUSY standard model are given in Eq. Q4.11| ). 
Thus, in our case, the Yukawa coupling constants of the lower generations will be renor- 
malized from the SUSY breaking to the GUT-scales by 

13/198 / ( s \ 3/2 / , x \ -8/9 



h u (M GUT ) 
f d (M GUT ) 
f e (M GUT ) 



h u (m z ) 
f d (rnz) 
f(%) 



^a 5 (M GUT ) 

ai{m z ) 
a b (M GUT ) 

a 5 (M GUT ) 



a 2 {m z 



a 5 (M, 



GUT 



7/198 / / \ \ 3/2 

/ a2{m z ) \ 



GUT) 



I a 3 {m z ) 
\a 5 (M GUT ) 

a 3 (m z ) \ 
a<j(M GUT )) 



-8/9 



,(B.6) 



3/22 / / \ \ 3/2 

| / OL 2 [m Z ) \ 



GUT) 



We have assumed the SUSY breaking scale to be m z . Here, h u , f d , f e are the Yukawa 
coupling constants of Hf and Hf to up-, down-type quark and charged-lepton multiplets. 
Note that the expressions differ from those given in Ref. The most important 

difference from the previous analysis is that the formula in Ref. does not have the 
factors of a 2 , and hence they overestimated the Yukawa coupling constants at the GUT- 
scale. 

Another ingredient which has not been included in Ref. is the large top-quark 
Yukawa coupling constant. Its effect appears in two points. First, the large Yukawa cou- 
pling constant contributes to the wave-function renormalization of the Higgs multiplet. 
If top quark becomes heavy, its large Yukawa coupling constant cannot be neglected in 
the RG equation of the Yukawa coupling constants of first and second generations. Its 
effect is to enhance the Yukawa coupling constants at the GUT-scale , and hence en- 
hancing the dimension-five operators. Second, the top-quark Yukawa coupling constant 

tThe authors in Ref. seem to have considered the renormalization of the mass operators rather 
than the Yukawa coupling constants. Since what determines the coefficient of the dimension-five op- 
erators is the Yukawa coupling constants themselves rather than the mass parameters induced by the 
SU(2) x U (1) breaking, it is clear that one should consider the running of the Yukawa coupling constants 
up to the GUT-scale. 
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itself become quite large at the GUT-scale, due to the wave-function renormalization 
of the top-quark multiplet. This further enhances the dimension-five operators for the 
third-generation contribution. Since the RG equation cannot be solved analytically if 
one includes the top-quark Yukawa coupling constant, we can only present the numerical 
results. 

From the Yukawa coupling constants at the GUT-scale, we obtain the coefficient 
of the dimension- five operators in Eq. ( |3.1| ). The dimension-five operators receive the 
anomalous dimensions, which we have to estimate to know their coefficients at the 
SUSY breaking scale. Here again, since the dimension-five operators are F-terms, all 
the renormalization effects come from the wave- function renormalizations of the external 
lines. We obtain the same result as that given in Ref. 0. The QQQL operators including 
only first- and second-generation fields in Eq. (|3.1|) receive an enhancement factor 

(B.7) 



, \ \ -1/33 / , \ \ —3 / , \ \ 4/3 

otxirnz) \ / a 2 {mz) \ / a 3 {m z ) \ 



\a 5 (M GUT ) J \a?,(M GUT )) \a 5 (M GUT ) J 

The enhancement factor of the operators including the third-generation fields is a little 
reduced by the effect of the top-quark Yukawa coupling. 

We show the numerical values of the coefficient As(i, j, k), the short-range renormal- 
ization effect between the GUT- and the SUSY breaking scales. First, if the contribution 
of the top-quark Yukawa coupling constant is neglected, this value becomes 



i \ \ 7 /" / ( \ \ ~ 4 / 9 



\a b (M GUT )) \a 5 (M GUT ) J 

0.59. (B.8) 



Thus our short range renormalization factor As is smaller than that in Ref. by | 
(their value is 0.91). Next, we show the numerical values of the coefficient As in Fig. 11 
for varying m t / a/2 sin (3h- Since As depend on the top-quark Yukawa coupling constant 
rather than m t , only a combination m t /y2 sin (3h is relevant. The solid line represents 
As for the dimension-five operators only with first- and second-generation fields, and the 
dash-dotted line for the operator {QtQt){QcL^}. The lower horizontal line is As with 
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the top-quark contribution neglected. One can see that the top-quark Yukawa coupling 
enhances the dimension-five operators. 

The factor Al in Eq. ( |3.12| ) is the long-range renormalization factor due to the QCD 
interaction between the SUSY scale and 1 GeV scale, and contain the renormalization 
of Yukawa coupling constants and the anomalous dimension of four fermi-operators. It 
is given in Ref. 0, 



A, 



/ q 3 (lGeV) r 2/3 f a 3 (m c ) \- 18/25 f a 3 (m b ) y 18/23 

0.22 (B.9) 



Combining all the renormalization effects, the four-fermi operators can be written down 
as in Eq. ( gl2|) . 



Appendix C Chiral Lagrangian Technique 

In this appendix we present a chiral Lagrangian technique for translating operators at 
the quark level to those at the hadron level. This technique has been developed in the 
Refs. 0, |2C| 

The baryon number violating four-fermi operators derived from the dimension-five 
operators are 



0{duuvi) 


= £ a p~f(dlu1)(dlv iL 


0(dudei) 


= e a ^(d a L ul){ule iL 


0(sudvi) 


= t a ^{s a L u p L ){dlv iL 


0(dusi , i) 


= e a fr{d a L ul)(slv iL 


O(suuei) 


= e aPl (slu L )(ule iL 



(C.l) 



where quark and lepton fields are in mass eigenstates, and i denotes the generation 
indices. There are also baryon number violating four-fermi operators derived from the 
dimension-six operators, and we will concentrate only on the operators relevant for the 
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decay mode p — > 7r°e" 1 , 

= e aM {dlui){ule L ), 

= e aPj (dt4)(ule R ). (C.2) 

The effective Lagrangian C q at the quark level for each decay modes can be written as 

C q {n,p — > ir{rf)i>i) = C(duub , i )0(duui' i ), 

C q (n,p — > Tr(r])ef) = C(dudei)0(dudei), 

C q (n,p — > Kui) = C \sudvi)0 (sudui) + C (dusu^O (dusui) , 

C q (n,p — > Kef) = C(suuei)0(suuei), 

C q {n,p^ne+) = CW6W + C®d<- 2 \ (C.3) 



The coefficients C's are derived from Eq. ( |3.12| ), and C's from Eq. ( A.l ) 



C(duuv,i) = — 2 . nQ 0- + V ){f(u, d) + f(u, e)), 

M Hc m l w sm 2(j H 

n<A a \ m u m d ^V: d y cd A L A s {u, c, dj) 

C(dud ei ) = — . (f(u,d) + f(d,p)), 

M Hc m w sm zp H 

C(sudui) = C{dusui) 

2ct\ m c m di e i ^V: d y cd V cs A L A s (c,u,d l ) tK 

= -T7 2 ■ no 0- + V )(f(u,d) + f{u,e)), 

M Hc 77% sm 2(j H 

n( \ 2a l rn u m d ^V: d y cs ALA s {u, c, dj) 

C(suue t ) = — , (f(u, d) + f(d, i/)), 

M Hc m w sm Zp H 



/7(i) _ p i^y5^R 

C-(2) = e ^4^(l + IK,l 2 



(C.4) 



We need to translate these effective Lagrangians at the quark level C q to the operators 
written in terms of the baryon and meson fields L h . 

Let us review the chiral Lagrangian for baryons and mesons. The Nambu-Goldstone 
bosons $ associated with the spontaneous breaking of chiral SU (3) l x SU (3)# symmetry 
can be written by 

U = expf^J (C.5) 
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where f w is the pion decay constant, and 



/ 



$ = 



TT 



K+ \ 



V 



TT 

K- 



(C.6) 



Similarly, the baryon fields can be written in the matrix form, 



B = 



( + V^ A 



P" 1 



n 



3 JV 



(C.7) 



Now the most general SU(3)l x SU(3)r invariant Lagrangian for strong interactions of 
mesons and baryons is 



\flTr(dU){drf) + Tr£(* ^ - M inv )B 



8 
+ 



+ l -i{D + F)Tr5 7 ^75 [C^C 1 ) - C^O] 5, 



where 



C = exp 



(C.8) 



(C.9) 



Since the (current) quark masses that break the chiral symmetry are small for up, down, 
and strange quarks, we can use Co to estimate the lifetimes of nucleon. 

Now we translate the effective Lagrangians containing quark fields C q to the ones 
containing baryons and mesons C h , by comparing the transformation properties under 
the SU(3)l x SU(3)r symmetry. The transformation properties of the baryon number 
violating operators given above are, 



0{duuvj),0{dudei) as (8,1), 
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0(sudi>i),0(suuei) as (8,1), 

O(sudi'i) as (8, 1), 

(1) as (3,3*), 

(2) as (3*, 3). (CIO) 

Thus, the four-fermi operators translating as (8,1) can be expressed in terms of the 
baryon and meson fields with a dimensionful constant /3, 

C h (n,p -> Tr{n)vi) = f3C(duup i )p dL Ti\P 1 CB L C} + h.c, 

C h (n,p^7i( V )et) = f3C{dude i )e dL Tr[P 2 (B L C]+h.c., 

C h (n,p -> Kiti) = i3C(sudu i )p dL Tr[P 3 CB L ^} + ^(dusz/^TrfP^C 1 ] + h.c, 

C h {n,p^Kef) = ^Cisuue^TriP^B^] + h.c, (C.ll) 

and the ones translating as (3,3*) or (3*, 3) with a dimensionful constant a, 



C h {n,p 



ire 



a 



CWe dL Tr[P 6 (B L (] + aC^e dR Tr[P 7 ^B R ^} + h.c, (C.12) 



with the coefficients C's and C"s defined in Eq. (|C.4|) . In the above formulae, Pi are 
projection operators, 



Pi 



V 



Pz 



( \ 



v 1 

( ^ 

0-10 

v 

( ^ 

-10 





Ft 



V 



^ ^ 



1 o o t 

^ 



o 1 1 

^ 



v 1 

( ^ 



1 



(C.13) 



/ 
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The constants a and (3 are the same as those defined in Eqs. (|A.3| ), ([3.19|) flj~9| . When 
one estimates the decay rates of the nucleons, one needs to include the virtual baryon 
exchanges. For example, to estimate the lifetimes of decay modes n,p — > Ki>, one 
should add the contributions from diagrams with virtual £ and A exchanges. We show 
the chiral Lagrangian factors of each decay modes in Table 5. We took ms = m-£ = 
m A = 1150 MeV, D = 0.81, F = 0.44 in Tables 1-3. 

Here we have some comments on the chiral Lagrangian factors in the decay rates of 
the nucleons. First, the ratio of Kv decay rates of neutron and proton is 

2 



r(p ->• K+V) 



2 _i_ 2m n f? 
m B 



1 + ^ {D + F) 



(C.14) 



This shows that the decay rate of the neutron is larger than that of the proton. The 
situation is different in ttu mode. The decay rate of the mode p — > 7r + z/ is two times 
larger than that of n — > 7r°z/. However, we use n — > 7r°z/ in section 5 because the 
experimental lower bound on r(n — > tt°u) is longer than that on r(p — > 7r + z/). 

Second, the decay rates into r\ of dimension-five operators are as large as those into 



7T 



For example, the ratio in v decay mode is 



r(n — > t)v) 

T(n — > 7T°Z/) 



m2) 2 3 



1 - UD-3F) 



\1 + D + F\ 



0.35. 



(C.15) 



Recall that the decay rates into rj from dimension-six operators is negligible fll9 |. Thus 
this mode may be interesting. 

The results in this appendix depend on the parameters a and /3. However, they are 
sensitive to hadron dynamics, and they differ for each hadron models. Even among the 



lattice calculations, the results vary from 0.03 GeV J IH31 to 0.0056 GeV J [E2L 
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Table 1 



The prediction of the nucleon partial lifetimes for the dominant decay modes, arising 
from the dimension-five operators (Q c Qc)(Q u L^) and {QtQt){QuL^). This class of decay 
modes depends on the parameter y tK . The mass degeneracy rric = mj and rrijx = m e - 
is assumed. The function / is defined in Eq. (|3.S| ) and depends on the SUSY particle 
masses. See the text for other variables. 

r(p -» K+Vn) = 6.9 x 10 31 
r(p -> K + v e ) = 1.4 x 10 33 

x 

r(n /f°P M ) = 3.9 x 10 31 
r(n -► i^°i? e ) = 7.7 x 10 32 

Table 2 

The prediction of the nucleon partial lifetimes for the next-leading decay modes, arising 
from the dimension- five operators {Q c Qc){QuL^) and {QtQt){QuL^). This class of modes 
depends on the parameter y t7T . The mass degeneracy = vn^ and = mg is assumed. 
The function / is defined in Eq. ( |3.8| ) and depends on the SUSY particle masses. See 
the text for other variables. 

r(p -> 7r+P M ) = 1.4 x 10 32 
r(p -> vr+Pg) = 2.9 x 10 33 

0.003 GeV 3 0.67 sin 2(3 H M Hc TeV" 1 2 

P A s 1 + y** 10 17 GeV f(u, d) + f(u, e) yrS 

r(n -> r/%) = 8.2 x 10 32 
r(n -> ?? z? e ) = 1.6 x 10 34 



0.003 GeV 3 0.67 sin 2(3 H M Hc 



TeV" 







A s 1 + y tK 10 17 GeV f(u, d) + /(«, e) 



yrs 



r(n — > tt°u, 



r(n — > TT^Uf 



2.9 x 10 



32 



5.7 x 10 



33 
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Table 3 



The prediction of the nucleoli partial lifetimes for the decay modes which depend neither 
on the parameter y tK nor y t7T . The relevant dimension-five operator is {Q u Qu){QcL 



alone. The mass degeneracy rric = vriu and mp = m^ e is assumed. The function / is 
defined in Eq. (|3.8|) and depends on the SUSY particle masses. See the text for other 
variables. 



r(p -» K°fi + ) = 1.0 x 10 35 
r(p -> 7r°/i+) = 2.0 x 10 35 



x 

r(p -> rj n + ) = 5.7 x 10 35 
r(n -> 7r-/i+) = 9.9 x 10 34 



0.003 GeV 3 0.67 _ s _ no M Hc TeV" 



■ sin 2(3h 



(3 A s 10 17 GeV f(u, d) + /(d, v) 



2 



yrs 



Table 4 

The experimental lower bounds on the nucleon partial lifetimes at the 90% C.L. WE 



r(p - 


♦ K+P) 


> 


1.0 x 10 32 yrs 


r{n - 




> 


8.6 x 10 31 yrs 


r(p - 


-* tt + D) 


> 


2.5 x 10 31 yrs 


r(n - 


-> 


> 


1.0 x 10 32 yrs 


r(n 


-> 77P) 


> 


5.4 x 10 31 yrs 


r(p -> 




> 


1.2 x 10 32 yrs 


r(p - 


> vrV + ) 


> 


2.7 x 10 32 yrs 


r(p - 




> 


6.9 x 10 31 yrs 


r(n — 


► 7T~/i + ) 


> 


1.0 x 10 32 yrs 


r(p - 


♦ vr°e + ) 


> 


5.5 x 10 32 yrs 
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Table 5 



Chiral Lagrangian factors in the nucleon-decay matrix elements. For notations, see the 
text. 



r(p - K + vi) 

r( n if ^) 
r(p tt+^) 

r(p - K°e+) 

r(p - vr e+) 

r(p ??e+) 
T(n - vr°e+) 
r(p 7r°e + ) 



(rrj - m|) 2 
32vrm3/2 

(m 2 - m^) 2 
32vrm3/2 



C(sudvi) — .D + Cidusvi 

6m b 



l + ^L(D + 3F) 
3m b 



C {sudfi 



1 



3ms 



(D - 3F) 



C{dusui 



m 



32vr/2 



^ \C(duuvi) [1 + D + Ftf 



m 



^-,\C{duuvi) [l + D + Ftf 



64tt/2 
(m 2 -m 2 ) 2 , 

647rm 3 j/ 2 
(m 2 - m|) 2 

32vrm3/2 



C(duuvi 
C(suuei) 



1--(D-3F) 



m. 



1-^E.m-F) 
m B 



m. 



'^\C{dudei)[l + D + F]Y 



64tt/ 2 
(m 2 - m 2 ) 2 
64vrm3/2 



C (dudei 



1--(D-3F) 



m 



327T/ 2 

m. 



^\C(dudei)[l + D + F]f 



64vr/, 



^(|c( 1 f + C( 2 ) 2 )[1 + J D + F] 



l + -JL(D + 3F) 
3ms 
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Figure Captions 



Fig. 1 A supergraph contributing to the dimension-five operators of the nucleon decay. 

Fig. 2 Allowed ranges on the color-triplet Higgs mass Mh c and the "GUT-scale" Mqut = 
(MyMs) 1 / 3 obtained from the renormalization group analysis (thick lines), by 
varying m~ h and m~ g between 100 GeV and 1 TeV. Mh c depends only on m~ h) and 
Mqut only on rrig. We use the gauge coupling constants at the weak-scale given in 
the text. Also shown are the ranges with an improved measurement on the strong 
coupling constant, a 3 = 0.118 ± 0.0035 (thin lines). 

Fig. 3 Lower bound on Mh c derived from the nucleon-decay experiments. The horizontal 
axis represents |l+y* x |, the sum of the second- and third-generation contributions 
normalized by the second-generation one. The vertical axis corresponds to M Hc . 
The shaded region is excluded. The upper curve corresponds to the hadron matrix 
element (3 = 0.03 GeV 3 , the lower one to (3 = 0.003 GeV 3 . The experimental limits 
come from the mode n — > K°P^ for |1 + y tK \ > 0.4, and from the mode n — > 7r°z/ M 
for |1 + y tK \ < 0.4 and |1 + y tw \ = 1. The short-range renormalization factor As 
is taken to be As = 0.67. The maximum value on Mh c {= 2.3 x 10 17 GeV) from 
the renormalization-group (RG) analysis requiring gauge coupling unification (see 
section 4) is also shown. 

Fig. 4 The dependence of the lower bound of Mu c on the parameters itLq and m^. The 
dashed line shows the dependence on taking uiq = 1 TeV. The dash-dotted line 
shows the dependence on tjiq when rriw = 45 GeV. In both curves we have taken 
the most conservative set of parameters, tan(3 H — 1, |1 + y tK \ < 0.4, |1 + y tw \ = 1, 
A s = 0.67, and (3 = 0.003 GeV 3 . We have assumed ~ itlq. The maximum 
value on M Hc {= 2.3 x 10 17 GeV) from the renormalization-group (RG) analysis 
requiring gauge coupling unification (see section 4) is also shown. 

Fig. 5 The dependence of the lower bound of M Hc on tan (3 H . The upper curve is obtained 
with |1 + y tK \ = 1, and the lower curve with |1 + y tK \ < 0.4 and |1 + y t7r \ = 1. We 
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have taken m Q = m L = 1 TeV, = 45 GeV, A s = 0.67 and (3 = 0.003 GeV 3 . 
The maximum value on Mu c (= 2.3 x 10 17 GeV) from the renormalization-group 
(RG) analysis requiring gauge coupling unification (see section 4) is also shown. 

Fig. 6 The limits on and rriq from the KAMIOKANDE nucleon-decay experiments, in 
the absence of the cancellation between second- and third-generation contributions 
(i.e., |1 + y \ = 1). The most conservative parameters, Mh c = 2.3 x 10 17 GeV, 
tan (3 H — 1, (3 — 0.003 GeV 3 , and A s = 0.67, are used. We have assumed nig ~ 
for simplicity. The shaded region is excluded. Also shown are the limits from the 
direct search experiments on wino and squarks at LEP and CDF. 

Fig. 7 The same as in Fig. 6, but allowing the cancellation between second- and third- 
generation contributions (i.e., |1 + y tK \ < 0.4, |1 + y t7T \ = 1). 

Fig. 8 The same as in Fig. 6, but with an improved constraint by a factor of 30 expected 
at superKAMIOKANDE. The expected limit on from the LEP-II experiment 
(rriu, > 90 GeV) is also shown. 

Fig. 9 The same as in Fig. 7, but with an improved constraint by a factor of 30 expected 
at superKAMIOKANDE. The expected limit on from the LEP-II experiment 
(rriyj > 90 GeV) is also shown. 

Fig. 10 The same as in Fig. 9. We have assumed that the error-bar of is reduced by 
a factor of 2 with the same central value, leading to a stronger upper bound on 
M Hc (< 6.1 x 10 16 GeV) (see Fig. 2). 

Fig. 11 The renormalization factor As vs. mt/\/2 sin/3//. The solid line represents the As 
for the dimension-five operators only with first- and second-generation fields, and 
dash-dotted line for the operator (QtQt)(Q c L^). The upper horizontal line is ^4^ 
derived by the authors in Ref. [0. The lower horizontal line is As which does not 
contain the contribution of the top-quark Yukawa coupling (i.e., m t = 0). 
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